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ON ESTIMATING GRAVITY ANOMALIES 
FROM GRADIOMETER DATA 


P. Argentiero 
R. Garza-Robles 

Goddard Space Flight Center 


INTRODUCTION 

The Gravsat-gradiometer mission as described in Reference 1 involves flying a gradioineter 
on a gravity satellite (Gravsat) which is in a low, polar, and circular orbit. The low-frequency 
terms of the geopotential field are observable from satellite perturbations as viewed by a 
high satellite or satellites. If the low-satellite orbit has an altitude of 300 km, a global dis- 
tribution of satellite-to-satellite tracking data can be obtained in about 10 days. From such 
a data set, an accurate and well-conditioned determination of the low-order terms of a 
spherical harmonic expansion of the geopotential field should be obtainable (Reference 2).* 
Because it measures second derivatives of the geopotential field, the gradiometer will be sen- 
sitive to the density variations of the Earth’s outer crust, which are the source of the high- 
frequency components of the geopotential field. These high-frequency terms tend to be 
smoothed out in satellite perturbation data, and, in this sense, the gradiometer as a data type 
complements rather than competes with satellite-to-satellite tracking data. 

The difficulty in using standard-estimation techniques to determine global geopotential fine 
structure is that a large number of parameters must be estimated. For instance, if features 
of the geopotential as small as 3 degrees are to be recovered and if the standard spherical 
harmonic expansion of the geopotential field is used as a parameterization, a full set of co- 
efficients to degree and order 60 is required. This implies the estimation of over 3700 para- 
meters. In general, simultaneous estimation of such large parameter sets is not possible. In 
practice, use of standard estimation techniques to recover geopotential fine structure from 
gradiometer data requires the adjustment of small subsets of parameters and constraining 
other parameters to a priori values. But, unless the parameterization exhibits a certain ortho- 
gonality property in the data set, the uncertainities in the unadjusted terms will badly corrupt 
the estimates of the adjusted terms. If spherical harmonic coefficients are used, the estima- 
tion problem can be decomposed into estimation problems of smaller dimensionality only 
if the gradiometer data set is globally distributed. This implies a heavy computational load. 

A parameterization of the geopotential field is desired in which the individual parameters 
exhibit a localized observability pattern in gradiometer data. With such a parameterization, 


*D. Koch and P. Argentiero, Simulation of the Gravsat iGeopaiise Mission, NASA TM X-70776. Submitted for presentation 
at the American Geophysical Union Fall Annual Meeting, December 12-17, 1974, San Francisco, California. 



it would be possible to accurately estimate small subsets of parameters by processing localized 
blocks of gradiometer data. The mean free-air gravity anomaly-Stokes’ function parameter- 
ization of geopotential fine structure (Reference 3) possesses a measure of orthogonality 
in local blocks of altimeter data (Reference 4). Reference 5 demonstrated the observability 
of gravity anomalies in gradiometer data, but the orthogonality property was assumed rather 
than demonstrated. Realistically, the orthogonality is not perfect and only the estimates of 
gravity anomalies which are a sufficient distance from unadjusted gravity anomalies will be 
of value. This implies that gravity anomalies should be estimated in blocks with the esti- 
mates of anomalies in a sufficiently small inner core accepted as valid and the rest should 
be rejected because of aliasing. To obtain an intelligent estimation strategy for determining 
gravity anomalies from gradiometer data, it is necessary to determine for a given data-block 
size the relationship between the accuracy with which a given anomaly is estimated and its 
distance from the nearest unadjusted anomaly. 


In this document, Gravsat-gradiometer missions at altitudes of 250 and 350 km are simu- 
lated, and optimal data-block and estimation-block sizes for determining 3- and 5-degree 
gravity anomalies are obtained by covariance analysis. A rotating-type gradiometer as 
described in Reference 6 is assumed to be mounted on the Gravsat and to have an accuracy 
of 0.1 etvos unit. However, the results can be readily scaled to reflect some other accuracy 
level. 


The rotating gradiometer is an instrument which senses gravity gradients in a given plane. 

If the instrument is onboard a satellite, that plane is assumed to be the orbital plane of the 
satellite. A local satellite-centered coordinate system (Ij , I^ , I3) is defined where the Ij 
unit vector is pointing northward, the I^ unit vector is pointing eastward, and the I3 unit 
vector is perpendicular to the plane spanned by Ij and I^ and is directed outward from the 
earth. If the satellite on which the rotating gradiometer is mounted is in a polar orbit, the 
sensing plane of the instrument always coincides with the plane spanned by the I^ and I3 
unit vectors. This assumption simplifies the following mathematical development. The 
geopotential field can be represented as 


W = U + T 


( 1 ) 


where U is a reference geopotential generally defined by a low degree and order spherical 
harmonic expansion of a nominal field and T is the so-called anomalous potential. The 
anomalous potential represents the high-frequency part of the field, and the primary 
interest is the portion of the gradiometer output which is attributable to T. Assuming 
that the satellite is in a polar orbit and neglecting off-diagonal elements of the gravity- 
gradient tensor which are several orders of magnitude smaller than the diagonal elements 
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(Reference 5), the additive component of the rotating gradiometer signal amplitude which 
is attributable to the anomalous potential T is 


AMP 


d2T(r,0, X) 1 a2T(r,0,A) 
9r^ 30^ 


( 2 ) 


where r, and X have the usual spherical coordinate interpretations and the derivatives 
are evaluated at the satellite position (Reference 5). 

The anomalous potential T at any point above the surface of the earth can be expressed 
by the discrete form of Stokes’ formula as 

T(r, 0. X) = 5gj(0',X')S(r. 0, X,0'.X')cos(0')A0[AX' (3) 

i 

where Sg. (0.', X.') is a mean gravity anomaly over a block centered at latitude 0! and longi- 
tude X! and referenced to a geoid defined by the nominal field U. The expression cos 
(0|) A0! AX! represents the area of the block on which the gravity anomaly is defined 
and S (r, 0, X, 0|, X|) is Stokes’ function given in Reference 3 by 


S(r, 0, X, 0', X') 



- 3 D - t cos ( 0) 


5 + 3 


1 + 1 cos ( 0) + D 


1 

}\ 


(4) 


where 


0 = cos"^ [sin (0) sin (0') + cos (0) cos (0') cos (X' - X)1 

^ ^ R 
r 


D = ( ] - 2 t cos (0) + )‘^“ 


where R is the mean radius of the Earth. The functional relationship between the anoma- 
lous part of the gradiometer output and a globally distributed set of gravity anomalies is 
obtained from equations 2 and 3 as 


AMP = 


3^S(r, 0, X, 0', X'-) I 3“S(r, 0, X, 0', X'.) 


ar^ 


302 


cos(0')A0'. AX;5g(0;,X') (5) 
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Equation 5 shows that, if the discrete form of Stokes’ formula is valid, the gradiometer 
output is a linear function of a global distribution of gravity anomalies. Therefore, stan- 
dard parameter-estimation techniques should apply to the problem of determining gravity 
anomalies from gradiometer data. To understand the application of these techniques, it 
is helpful to use matrix notation. Let y be a vector of the anomalous components of 
gradiometer readings. Assume that the i^^ component of y is the anomalous gradiometer 
reading obtained when position of the satellite is given in spherical coordinates as (r., 0., 
X.). Let g^be a vector of numerical values of a global set of gravity anomalies. The 
component of g" is the numerical value of the gravity anomaly centered at latitude 0f and 
longitude X!. The functional relationship between y” and g" is expressed as 

'y = Ag' (6) 

where A is a matrix the number of whose rows is the number of observations and the 
number of whose columns is the number of gravity anomalies. The element in the i^^ row 
and column of A is 


A(l,j) = 


a^Slr, 0i,Xi,0j,Xj) 


1 a 


•S(ri,0i.Xi,0'.Xp 


3rf 


90r 


cos (0j) A0j AXj (7) 


Equation 7 provides a linear equation of condition and, in a standard minimum variance 
fashion "g, could be estimated from observations of y. In order for equation 7 to be correct 
(correct, that is, assuming that the approximations inherent in the discrete form of Stokes’ 
formula are valid), the gravity anomalies in the array g"must cover the globe. Computa- 
tional considerations make it necessary to choose a region in which the number of gravity 
anomalies to be estimated does not exceed a few hundred. In effect, gravity anomalies 
outside of this region are assumed to be zero. To see precisely what happens when this 
assumption is made, postulate that the g" array of equation 6 is to be defined over the globe, 
and write 


g = 





h 


where 

g^j = gravity anomalies to be adjusted in a standard minimum-variance filter and 

'g^ = gravity anomalies assumed to be zero and therefore left unadjusted by the 
minimum-variance filter 

Then equation 6 can be written 


y = A|g, + Ajgj 


( 8 ) 
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where Aj and are the variational matricies of y with respect to and , respectively. 
The gradiometer output provides direct observations y of y with statistics 

y = 'y + E (i^) = 0, E = Q (9) 

Estimates of mean free-air gravity anomalies obtained from satellite tracking and gravi- 
metry measurements are available (Reference 7). Unless this information was correctly 
factored into the gravity anomaly estimates obtained from gradiometer data, the resultant 
estimates would not be optimal. Consequently, the existence of an a priori estimate g' or 
is assumed with statistics 

g' = 'g J + ocj , E (aj ) = 0, E (a^ o:[ ) = (10) 

For computational reasons, the gravity anomalies g^ are assumed to be zero, but the 
actual values of gravity anomalies in the region of the sphere which is ignored have a 
certain distribution about zero. Assume that 

E (g^ ) = 0, E (g, g] ) = Pj (11) 

When the values of g^ are assumed to be zero, the minimum variance estimate of^j 
becomes 

g, = (A[Q-'A| +P-')-' (A|Q-'y+ P-Ig',) (12) 

Define the covariance matrix of the estimator given by equation 12 as 

P = E(|2| -g| 1 Ig, 1 (13) 

From equations 9, 10, 11, and 12, 

g, -g, = (A|Q->A, +P-')-> (-A|Q-' A,g2 + AjO-'i' + p-'a,) (14) 

Equation 14 yields 

p = (A[q-' a, + p-' )■' +(AJ q-' a, + p;' )■' A[Q-' A,P,AJQ"' A,(Af Q-' A, + P;' )-' (15) 

Assume that the data noise is uncorrelated and that each data point has the same variance. 
Then 

Q = I<^o C6) 
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where I is the identity matrix and is the common variance of the data. On the assump- 
tion that the values of the unadjusted gravity anomalies are independently distributed, the 
covariance matrix of can be written as 


■> 



^2 = 


0 


(17) 



where n^ is the number of unadjusted gravity anomalies and a? is the second moment about 
zero of the i^^ unadjusted gravity anomaly. Also, define a matrix K as 

K = (A[Q-' A, +P-‘ )■' A[Q-'A2 (18) 

If iij is the number of adjusted parameters, then K is of dimension iij by . With these 
assumptions, equation 1 5 yields the following expression for the variance of the i''* adjusted 
gravity anomaly 


2 

P(I, I) = j aj)2 (19) 

j = 0 

where /3? ^ is the i*''‘ diagonal element of the matrix (Aj^ Aj )"' (assuming that diagonal 
elements of the matrix Pj^ are relatively small) and 

= K(I.j)J> 1 (20) 

Define the error sensitivity matrix as 

S - K2 n,.J = 0, 1,2 n, (21) 

And finally define the alias matrix as 

L = Sa (22) 
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where 


o 


0 


o = 


0 



a 


n 


2 


(23) 


The alias matrix reveals much of the probability structure of the estimation procedure. 
From equations 19, 21 , 22, and 23 it can be seen that the standard deviation of the 
adjusted gravity anomaly is the root sum square (RSS) of the terms in the i^'^ row of the 
alias matrix. The elements in the first column of the alias matrix represent the RSS contri- 
bution to the standard deviation of each estimated parameter due to the data noise. The 
elements in the column,] > 2 represent the RSS contribution to the standard deviation 
of each estimated parameter attributable to the] - 1®^ unadjusted parameter. These terms 
are called the aliasing contributions to the uncertainty in the adjusted parameters attribut- 
able to the uncertainty of the] - unadjusted parameter. Note that the aliasing contri- 
butions attributable to the j^'^ unadjusted parameter are proportional to the standard 
deviation of the j^* parameter. 


DEFINITION 

In a given estimation process, the i^'* estimated parameter is 
said to be orthogonal with respect to the j^^ unestimated 
parameter if the aliasing contribution to the i^'* estimated 
parameter attributable to the uncertainty of the j^'* unesti- 
mated parameter is zero. It will be seen that the relation is 
symmetric in the sense that, if one parameter is ortliogonal 
with respect to another and if they are interchanged within 
the adjusted and unadjusted modes, they are again orthog- 
onal. 

To understand the implications of the orthogonality relationship, a more revealing repre- 
sentation of the aliasing terms is necessary. Note that the first term on the right side of 
equation 1 5 is the covariance matrix of the estimation process under the assumption that 
the unadjusted parameters are perfectly known. This covariance matrix gives the uncer- 
tainty of the estimates which are attributable to the data noise only. Define the so-called 
“noise-only” covariance matrix as 

P = (A| Q-' A, + P-i )-' (24) 
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Next, observe that the elements in the row and column of and , respectively, 

are the partial derivatives of the i''’ data point with respect to the adjusted parameter 

and the partial derivatives of the data point with respect to the unadjusted parameter. 

The aliasing contribution to the adjusted parameter attributable to the j^'* unadjusted 
parameter can be written as 


L(IJ + 1) 


m 


E E 

k = 1 e= 1 


ay(g) 

ag,(k) 


Q-'(C,C) 


3y(g) 

0g2(j^ 


(25) 


where m is the number of data points. If the estimates of the adjusted parameters are 
relatively uncorrelated in the noise-only covariance matrix, equation 25 can be approxi- 
mated by 


m 

Ldj + 1) = P(i, 1) ^ 
t = 1 


9y(g) 

ag,(i) 


Q-'(e. 0 


dyie) 




(26) 


Equation 26 shows that the aliasing contribution is symmetric in its arguments, and this 
implies the symmetry of the orthogonality relationship. 

A sufficient condition for the left side of equation 25 to approximate zero is for the 
observability patterns of g^ (I) and (j) in the data to be virtually nonoverlapping. Fig- 
ures 1 and 2 demonstrate that the perturbation of a gradiometer observation caused by an 
0.01-mm/s^ (1-mgal)* perturbation of a gravity anomaly rapidly attenuates with increasing 
spherical radius. Therefore, if the grids on which (I) and g^ G) are defined are sufficiently 
separated, the orthogonality relationship will be effectively satisfied, and the estimate of g'^ (I) 
would experience no aliasing from the uncertainty of g^G)- Conversely, if the grid on 
which gj (I) was defined was close to grids whose gravity anomalies were adjusted, serious 
aliasing of the resultant estimate would probably occur. 

Therefore, if the gravity anomalies are estimated in a block, the outer layers of the block 
contain gravity anomalies whose estimates will be badly aliased by the adjacent unadjusted 
parameters. It will be necessary to discard these estimates. However, the gravity anomalies 
in a sufficiently small inner core of the block may be adequately separated from the unad- 
justed parameters to be effectively orthogonal with respect to them. The estimates of 
these terms will presumably be of sufficient accuracy to be accepted. In effect, for every 
block of gravity anomalies to be estimated, it will be necessary to construct a “buffer 
zone” several layers deep of gravity anomalies which surround the block. The new and 


*TIiroughout the text of this document, the measurement unit of miUigalileo (mgal) has been converted to a Standard Inter- 
national Unit of miUimeters/sec^ (mm/s^). For convenience, the illustrations have not been converted. The simple con- 
version is 1 mgal = 0.01 mm/s^. 
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SPHERICAL DISTANCE 

Figure 1. Gradiometer output perturbation attributable to 
0.01 -mm/s^ (1-mgal) perturbation of 3-degree gravity anomaly. 

larger block of gravity anomalies must be simultaneously estimated and because of the 
estimates of gravity anomalies in the buffer zone must then be rejected aliasing. To design 
an intelligent data-processing procedure for the Gravsat-gradiometer mission, the relation- 
ship between the depth of the .buffer zone and the accuracy of the estimation procedure 
must be determined. This relationship will vary with grid size, data-block size, and satellite 
altitude. The most efficient way to study the relationship is to use, covariance-analysis 
techniques to generate alias matricies for several situations and to attempt generalizations 
from the results. This method is explained in the next section. 

RESULTS 

Assuming that Gravsat altitude is approximately 300 km and that a gradiometer reading 
is obtained every 15 seconds, a dense, globally-distributed set of approximately 56000 
gradiometer measurements is obtained within 10 days. In the low latitudes, the subearth 
points of the satellite at the measurement times form a grid on the Earth of 1 degree lati- 
tude by 1-1/2 degrees longitude. All simulations rely on this assumption concerning the 
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SPHERICAL DISTANCE 


Figure 2. Gradiometer output perturbation attributable to 
0.01-mm/s^ (1-mgal) perturbation of 5-degree gravity anomaly. 


data distribution. The dynamical aspect of the problem is therefore eliminated, and the 
aliasing effect of orbit-determination errors are not included in the results. This shortcoming 
is not serious. A Jet Propulsion Laboratories (JPL) study (Reference 8) indicates that, 
for an instrument with an accuracy of 0.1 etvos unit, orbital accuracies of 20 to 25 meters 
are adequate. Satellite-to-satellite tracking of Gravsat from a high relay satellite should 
yield orbit-determination accuracies well within this range (References 2 and 9). 

.As mentioned in the introduction, the reason for using gravity anomalies to describe geo- 
potential fine structure is the possibility of successfully estimating local blocks of gravity 
anomalies in local blocks of gradiometer data. To investigate this possibility, attention 
is focused on the estimation accuracy of a single gravity anomaly as a function of the 
number of adjacent gravity anomalies which are simultaneously adjusted and the quantity 
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of data which is used in the estimation. The estimation block is defined as a block of 
gravity anomalies centered at the gravity anomaly in question and square in the sense that 
the block subtends the same number of degrees longitude as latitude. All the gravity 
anomalies in the estimation block are assumed to be simultaneously adjusted in a given 
estimation procedure. As defined, estimation blocks can be incremented only by integral 
numbers of twice the gravity-anomaly block size. For 3-degree gravity anomalies, for 
instance, possible estimation block sizes are 3, 9, 15,21 degrees, etc. 

The data block is defined as a gradiometer data set which was obtained when the sub-Earth 
points of the satellite were within a square region of the Earth’s surface centered at the 
gravity anomaly in question. For a given estimation-block size, the accuracy of a gravity 
anomaly estimate is not necessarily a monotone function of the data-block size. At first, 
this may be puzzling because better estimates are generally expected from increased data. 

An explanation can be obtained by noting that for a given estimation block the covariance 
matrix of the set of estimated gravity anomalies is given by equation 1 5 as the sum of a 
matrix which is dependent only on data uncertainty and a matrix which represents the 
aliasing effects from the unadjusted parameters. With increasing data blocks, the elements 
of the first matrix must decrease, but the elements of the matrix which conveys the alias- 
ing effects will generally increase. This effect can be shown graphically by means of so- 
called aliasing maps. To obtain the aliasing maps of figures 3 and 4, in a covariance mode 
a situation was simulated in which the geopotential fine structure was described by means 
of 3- by 3-degree gravity anomalies and Gravsat-gradiometer data with an accuracy of 
0.1 etvos unit was assumed. In figure 3, the RSS contribution was mapped in mgals to the 
uncertainty in the estimate of the gravity anomaly defined on the blank grid element in 
the middle of the graph when the adjacent gravity anomalies were assumed to be in an 
unadjusted mode and to have an a priori uncertainty of 0.50 mm/s^ (50 mgal). The data 
block represented by the shaded area of the graph was 10 degrees on a side. The satellite 
altitude was 250 km. Note that the aliasing contributions decrease with the distance between 
an unadjusted parameter and the estimated parameter, thus demonstrating the inherent 
orthogonality properties of gravity anomalies in gradiometer data. In figure 4, the data- 
block size has been increased to 20 degrees, and the aliasing effect of adjacent unadjusted 
gravity anomalies is greater. 

The choice of an intelligent data-reduction strategy is obviously dependent on a knowl- 
edge of the relationship between estimation accuracy and the choice of data-block size 
and estimation-block size. A computer program was written which computes the right 
side of equation 15 for any given data-reduction technique. The program was used to 
investigate data-reduction strategies. First, an altitude for Gravsat was chosen and a gradi- 
ometer data block centered on a chosen gravity anomaly was assumed. The data accuracy 
was assumed to be 0.1 etvos unit. A normal matrix for 400 gravity anomalies was formed, 
Tlie chosen gravity anomaly and gravity anomalies in successive layers surrounding the 
chosen gravity anomaly were assumed to be in an adjusted mode, and the rest were placed 
in an unadjusted mode with uncertainties about zero of 0.50 mm/s^ (50 mgal). The right 
side of equation 15 was computed, and the resultant standard deviation in the estimate of 


11 



10° DATA BLOCK 
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Figure 3. Alias map for 3-degree gravity anomalies. 

the chosen gravity anomaly was recorded. The estimation-block size was incremented by 
removing another layer of gravity anomalies from the unadjusted to the adjusted mode and 
repeating the process. Figures 5 and 6 show the relationship between estimates of 3- and 
5-degree gravity anomalies and the choices of data-block size and estimation-block size when 
Gravsat is at an altitude of 250 km. Figures 7 and 8 show the same information when 
Gravsat is at 350 km. The graphs show that an intelligent compromise between computa- 
tional load and accuracy is achieved when a data-block size of 15 degrees is used and the 
buffer zone of rejected gravity-anomaly estimates is at least three layers deep. The graphs 
show that when this strategy is followed and when the satellite is in a 250-km orbit, 3- and 
5-degree gravity anomalies can be estimated with uncertainties of 0.03 and 0.01 mm/s^ (3 and 
1 mgal), respectively. When the satellite is in a 350-km orbit, 3- and 5-degree gravity anomalies 
can be estimated with uncertainties of approximately 0.07 and 0.025 mm/s^ (7 and 2.5 mgal), 
respectively. 

These results are dependent on the chosen value for the accuracy of the rotating gradiom- 
eter. Because this choice is somewhat arbitrary, it is useful to parametrically represent 
gravity-anomaly-estimation accuracy as a function of gradiometer accuracy. Equation 19 
implies that, for a given choice of estimation-block size and data-block size, the standard 
deviation of a gravity-anomaly estimate is the root sum square of a term which represents 
the aliasing effects of unadjusted gravity anomalies and a term which is proportional to 
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Figure 4. Alias map for 3-degree gravity anomalies. 

the standard deviation of the rotating gradiometer observation. Assuming a IS-degree data- 
block size and assuming that a given gravity anomaly is simultaneously estimated with all 
gravity anomalies in the three adjacent layers, the estimation accuracies for 3- and 5-degree 
gravity anomalies at satellite altitudes of 250 and 350 km are given as a function of data 
accuracies as 


a (3°. 250km) = (625 + 2.25)’'^ 

a (5°, 250km) = (81 + 0.5(10)--)’''^ 

a (3°, 350km) = (4200 a’ + 3.24)'/^ 
a (5°, 350km) = (625 + 0.01 )'/^ 


(27) 


Standard deviations on the left side of equation 27 are in mgal and the data standard devia- 
tion a on the right side of equation 27 is in etvos units. 
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Figure 5. Accuracy of 3- by 3-degree mean free-air gravity anomaly estimate versus data- 
block size for various estimation-block sizes. 



Figure 6. Accuracy of 5- by 5-degree mean free-air gravity anomaly estimated versus data- 
block size for various estimation-block sizes. 
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'igure 7. Accuracy of 3- by 3-degree mean free-air gravity anomaly estimate versus data- 
block size for various estimation-block sizes. 
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Figure 8. Accuracy of 5- by 5-degree mean free-air gravity anomaly estimate versus data- 
block size for various estimation-block sizes. 
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CONCLUDING REMARKS 


If a rotating gradiometer were flown on a satellite in a polar, circular orbit with an altitude 
of 250 to 350 km, a dense and globally distributed set of gradiometer observations would 
be obtained in about 10 days. From this data set, it should be possible to determine a global 
set of gravity anomalies. However, to keep the computational load within reasonable bounds, 
it will be necessary to estimate local blocks of gravity anomalies in local blocks of gradiom- 
eter data. The results of this study indicate that good gravity-anomaly estimates can be 
obtained from gradiometer data blocks as small as 15 by 15 degrees if each gravity anomaly 
whose estimate is to be accepted is separated from unadjusted gravity anomalies by three 
layers of adjusted gravity anomalies. If the satellite is flown at an altitude of 250 km and 
if the instrument has an accuracy of 0.1 etvos unit, this estimation strategy yields estimates 
of 3- and 5-degree gravity anomalies accurate to within 0.03 and 0.01 mm/s^ (3 and 1 mgal), 
respectively. When the satellite is flown at an altitude of 350 km, the accuracies of 3- and 
5-degree gravity-anomaly estimates are 0.07 and 0.025 mm/s^ (7 and 2.5 mgal), respectively. 
Intermediate results can be expected at intermediate altitudes. Equation 27 of this document 
permits a scaling of these results to a data accuracy level different from 0.1 etvos unit. 


Goddard Space Flight Center 
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Greenbelt, Maryland June 7, 1976 
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